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We calculate ground and excited electron and hole levels in spherical Si quantum dots inside 
SiC>2 in a multiband effective mass approximation. Luttinger Hamiltonian is used for holes and the 
strong anisotropy of the conduction electron effective mass in Si is taken into account. As boundary 
conditions for electron and hole wave functions we use continuity of the wave functions and the 
velocity density at the boundary of the quantum dots. 
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I. INTRODUCTION 

The study of materials composed of Si nanocrystals 
dispersed in SiC>2 matrix is an issue of high importance 
for various optoelectronic applications 1,2 . In particular, 
the knowledge of the energy spectrum of carriers confined 
in the nanocrystals and their wave functions is crucial 
for the understanding of electronic processes. Numer- 
ous theoretical works have been developed for the eval- 
uations of the ground-state electron-hole pair energy of 
Si nanocrystals^^i^. On the other hand, the prob- 
lem concerning the energy-level positions and the corre- 
sponding eigenfunctions of the excited carrier states is 
much less studied^^. However, shortly after the gener- 
ation of an electron-hole pair in a nanocrystal a number 
of important non-equilibrium processes involving these 
"hot" carriers may take placed necessitating thus the 
knowledge of the excited states. 

This current paper is devoted to the study of the 
ground and the excited electron and hole levels within 
the framework of a multiband effective mass approxima- 
tion. The finite energy barriers at the Si/Si02 bound- 
ary are explicitly accounted for. For the description of 
the electron and hole states of the carriers confined in Si 
nanocrystals we utilize the envelope function approxima- 
tion taking into account the elliptic symmetry of the bot- 
tom of the conduction band and the complex structure of 
the top of the valence band in Si. The finite energy bar- 
riers at the boundary between Si and SiC>2 are treated by 
employing the Bastard boundary condition a 12 i 13 . This is 
mainly the difference of our method comparing to earlier 
calculations based on the effective mass approximation 9 
that failed to describe properly the optical properties of 
small nanocrystals. An advantage of our method in com- 
parison with ab-initio methods based on the density func- 
tional theory^ is that without a considerable numerical 
effort we can calculate not only the ground state but also 
excited states of the confined carriers. Furthermore, our 
theory can be applied to a broad range of nanocrystal 
sizes. We are not limited to nanocrystals with a small 



number of atoms. The other point is that the calcu- 
lations of various excitation and de-excitation processes 
under participation of the confined carriers using the de- 
rived wave functions are transparent and allow for an 
insight into the underlying physics. 



II. ELECTRON STATES 

The conduction band of bulk Si has six equivalent 
minima in the first Brillouin zone at points iko. z = 
(G,0,±0.85)ftjc, ±k ,y = (0,±0.85,0)fc x , and ±k^ x = 
(±0.85, 0, 0)fcjtj where kx — 2n/a and a = 0.543 nm is 
the lattice constant of Sii£. The minima are situated in 
the neighborhood of the six A-points (there are three 
non-equivalent A-points). The conduction band is dou- 
bly degenerate at each of the A-points, which is a con- 
sequence of the fact that Si lattice has two atoms in the 
elementary unit cell and the origin can be chosen at the 
center of any of them. Assuming the Bloch amplitudes 
not changing in the neighborhood of the A-point one 
can write the wave function of one of the six equivalent 
ground states of electrons in the nanocrystal as 

C = f (r> c ,e^- f (y = ±x, ±y, ±z), (1) 

where u cu is one of two Bloch amplitudes of bulk electron 
at A-point in the Brillouin zone, which corresponds to 
the lower conduction band at the k$ tV point. Note, that 
u cv yields zero overlap integral with the Bloch amplitudes 
of the top of the valence bancU^. The overlap integral of 
the Bloch amplitudes of the top of the valence band with 
the second Bloch amplitude u 8V — u cl ,e~ %2kx -" r corre- 
sponding to the upper conduction band at ko tll is not 
equal to zero. The envelope wave function £ e in Eq. ((T|) 
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inside the Si quantum dot satisfies the following equation: 



h 2 d 2 

2m\\ dz 2 



2m i 



d 2 d 2 \ 
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+EC(x,y,z) = 0, 



(2) 



where mil = 0.916too, m± = 0.19too with too being 
the free electron mass. The rigorous formulation of 
the boundary conditions for the boundary between Si 
and SiC>2 in the framework of the envelope function 
method is not a trivial task and generally has to be in- 
vestigated in comparison with experiment and numer- 
ical methods. The Bastard type boundary conditions 
imply that £ and vt; are continuous across the bound- 
ary, where v = i[r, -H] is the velocity operator, £ is 
the envelope wave function and H is the correspond- 
ing Hamiltonia n 12 ' 13 . Here we assume that the spectrum 
of electronic states outside the nanocrystal in SiC>2 is 
isotropic and determined by a single electron effective 
mass, which is equal to the free electron mass mo- Then 
outside the quantum dot we have 



n 2 + {p_ + &_\ 

2too \dx 2 dy 2 dz 2 J ' ' 

+(E-U e )e(x,y,z) = 0, 



(3) 



where U e is the energy barrier for electrons. According 
to Ref.Qi^we have U e = 3.2 eV. The boundary conditions 
result in the following equations 
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FIG. 1: (Color online) Dependence of positions of the electron 
energy levels above the bottom of the conduction band of 
bulk-Si on the quantum dot diameter. The numbers near 
the lines indicate the total degeneracy (including the spin 
degeneracy) of the corresponding levels. 



the lower energies. The reason why we have smaller en- 
ergies is that this tight-binding model used truncation of 
Si nanocrystals by H atoms. This procedure is known to 
give higher energies and greatly overestimate the optical 
band gap when compared with experiments on Si/Si02 
nanocrystals and recent ab-initio TD-DFT calculations^ 
as well as with our model. 



III. HOLE STATES 



m± R n 



dp 



dz 
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- toq dr 



(5) 

which are most conveniently written in the cylindrical 
coordinate system (z, p, (j)). 

Equations ([2]) and ([3]) with the boundary conditions ((4)) 
and ([5]) have been solved numerically after separating the 
trivial angular part -y== exp(imcf>) (m = 0, ±1, ±2, ...) of 

the wave functions. We have got the electron energy lev- 
els and the corresponding envelope wave functions. The 
dependence of positions of the several lowest energy levels 
on the quantum dot diameter is depicted in Fig. [TJ The 
dependence of the electron envelope wave functions on 
the distance from the center of the quantum dot is shown 
in Fig.[5]for d — 2 nm (such a small diameter is chosen for 
demonstration reasons in order to resolve better the tun- 
nelling tails of the envelope wave functions). We have 
also compared the positions of the electron levels and 
their degeneracies with the existing data of Ref. calcu- 
lated by the empirical tight-binding method for quantum 
dots with diameter d — 7.61 nm. We have found that 
apart from small level splittings due to the valley-orbit 
interaction neglected in our model we get the same se- 
quence of levels. The levels are, however, shifted towards 



For description of the valence band structure in Si we 
use a generalization of the Luttinger Hamiltonian^ in the 
limit of a vanishing spin-orbit coupling, which is justified 
for Si, i.e. we write the Hamiltonian H in the form 



H = (A + 2B)p 2 - 3B(p ■ J) 2 



(6) 



where p is the momentum operator and J is the unitary 
angular momentum operator acting in the space of Bloch 
amplitudes. Furthermore, we introduced 



4 TO/jTO; 



mi 1 m h 
— , a = — 

4 m h mi 



rn, 



(7) 



m 



7i - 2 7 



TO; 



to 



7i + 2 7 



7 = 7( 3 73 + 272). (8) 



Values of the constants 71, 72, and 73 for Si are 4.22, 
0.53, and 1.38, respectively^. The basis of the Bloch 
amplitudes space can be chosen in the form of spherical 
components^ u a = Z, u± = =F^/l/2(X ± iY) of the 
corresponding functions X = yz, Y = xz, and Z = xy, 
of the representation T25— 

In the bulk Si this model leads to two types of states 
corresponding to a doubly degenerate (in absence of spin- 
dependent interactions) heavy hole band having mass 
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FIG. 2: (Color online) Electron wave functions in dependence on the cylindrical coordinates z (a) and p (b) for five lowest 
electron levels of a quantum dot with diameter of 2 nm. The wave function £J of the fourth level with magnetic number m = ±1 
has also an angular dependence e lmi *. 



rrih and a non-degenerate light hole band having mass 
2mhmi/(3mh — mi). The quantum confinement gives 
rise to mixing of the states. Eigenfunctions of the 
Hamiltonian ([6]) can be found as eigenfunctions ipFM of 
the square F 2 of the full angular momentum operator 
F = L + J [L = —irx dp) and its projection F z onto the 
axis Eigenvalues of F 2 and F z are F(F +1) and M, 
respectively, where F = 0, 1, 2, ... and M can be any inte- 
ger number having absolute value not larger than F. For 
a spherical quantum dot there are three types of states, 
namely 

KM' F+1 (r,e,q>) =R F F -\r)Y F F ll(6,cp) 

+ i?£+V)YSM0), 

V^ M (r, 8, 0) = R F (r)Y F M (9, </>), (10) 

^o(r,M) = i$(r)Yi,(M), (11) 

where i?£ -1 (r), i?^ +1 (r), and i?£(r) are the radial parts 
of the envelope wave functions. Furthermore, 

Yfw(M)= E ^l^WMK* (12) 
mi ,m2 

i 



are the vector spherical harmonics that can be expressed 
in terms of the usual spherical harmonics Y nm {6, (f>) and 
the Clebsh-Gordon coefficients C^"L „■ „, — . For the first 
two types of functions © and (flU|) only solutions with 
F ^ 1 are possible. We will see that the functions © 
are of a mixed type, whereas the functions (flTj|) and (TlTj) 
are of a heavy and light hole type, respectively. 

For a formulation of the boundary conditions for the 
hole states we take into account that the main contri- 
bution to the valence band states in SiC>2 is given by 
p-orbitals 20 and for description of hole states outside the 
Si nanocrystal we can use the same form of the Luttingcr 
Hamiltonian © . It is known that the hole masses at the 
valence band maximum in Si02 are pretty large^S. For 
simplicity we choose them to be equal to m v = 5too- 
Then the corresponding values of the coefficients A a and 
B Q of the Luttinger Hamiltonian © are A a = — and 
B Q = 0. In such a case we can formulate appropriate 
Bastard type boundary conditions. 

Inserting the functions V'fm 1 ' +1 given by Eq. ([9]) into 
the Schrodinger equation with the Hamiltonian we get 
the following equation system (see also Ref. KL8I ) for the 
radial functions i?^ _1 (r), R F+1 (r) inside the nanocrystal 
(r < R nc ): 
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where E denote the hole energy and fi = B/A. The 
general solution of the equation system of Eqs. (flU)) and 
(fT4|) . which does not diverge at r = 0, is found as 



Rf' 1 ^) = Cj F ^{Xr/R nc ) + Djp-iiXpr/R,*), (15) 
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(16) 



where C, T> are coefficients to be found from the bound- 
ary and normalization conditions, ji(z) are the spherical 
Bessel functions of the first kind 21 , and 



1 + 2// 



(17) 



The energy E, which is negative, is connected with the 
positive variable A via 



Ah 2 , . o 
£ = — (1 -^A 2 = - 



Rl c K ^ 2m h R 2 n 
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Outside of the nanocrystal (r > -Rnc) the radial parts 
ns 

12 2d {F-1)F~ 



of the functions "4> F m' F+1 satisfy the following equations 



A h z 



_d_ 

dr 2 



r dr 



(19) 



= -{E + U h )R F F -\r) 
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A ft 2 
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(20) 



-{E + U h )R F F + \r) , 



where f/ft is the energy barrier for holes at the Si/SiC>2 
boundary. According to Ref. [if] we have Uh = 4.3 eV. 
The general solution of the equation system of Eqs. (fl"9|) 
and (|20|) , which converges to zero for large distances from 
the nanocrystal is found as 
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where 



£ _1 (r) = (C + 2? )fc F _i( K r/i? nc ), 



(21) 
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(22) 



K = 



v / 2m v (E + U h ) R 



(23) 



C , T> are again coefficients to be found from the bound- 
ary and normalization conditions, and ki{z) are the mod- 
ified spherical Bessel functions of the third kind 2 ^ . The 
boundary conditions lead to the following equations for 
the radial functions: 



R F F -\r) 
R F F +1 (r) 



r=Rnc 



R F F - l {r) 
R F F +1 (r) 



(24) 
(25) 
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A f r R F +i (r) 



(26) 



r=_R+ 

(27) 



Using the functions (T5]), (fH) ]) and ([2"T] ) . ([22] ) in Eqs. 
([27]) leads to a solvability condition determining the en- 
ergy eigenvalues for states ip F M~' F+1 - We have derived 
this condition (see Appendix). From it we have found 
numerically the energy eigenvalues in dependence on the 
nanocrystal radius R nc . The corresponding coefficients 
C,T>, C , and T> , assuring the normalization condition, 
have been also derived numerically. 



For the radial function R F (r) of the states tp FM we get 
the following equation inside the Si quantum dot 



dr 2 



2 d 

r dr 



F(F + 1) 



E 



R F (r) = -^4(r). 

(28) 

One can easily see that it is the same equation as for the 
radial part of wave function of a particle having a simple 
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FIG. 3: (Color online) Dependence of the positions of the 
hole energy levels below the top of the valence band of bulk- 
Si on the quantum dot diameter. The numbers near the lines 
indicate the total degeneracy (including the spin degeneracy) 
of the corresponding levels. Types of the levels and values of 
the total angular momentum F are also shown. 



parabolic band with the heavy hole mass and angular 
momentum F. For the radial function i?J(r) of the states 

V'oo we g et 



(1 + 2*0 



dr 2 



2 d 

r dr 



RUr) = -jjpR&r) . (29) 



This equation is then the same as for a simple particle 
having the light hole mass and angular momentum 1. 
Solving equations f2"5]) and (|2"9"j) we find 

(30) 



R%{r)=C'j F (\r/R nc ), 



R 1 (r)=C"j 1 (\pr/R ac ), (31) 

for r < R nc , where C and C" are the corresponding nor- 
malization coefficients. Outside the nanocrystal Rp(r) 
satisfies the equation 
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and i?o(r) satisfies the equation 
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Rl(r) = -(E+U h )Rl(r). (33) 



Solutions having an appropriate behavior at infinity are 

B%{r)=C' o k F {Kr/R 0D ), (34) 

Rl{r)=C^k 1 {Kr/R ac ). (35) 
The boundary conditions in this case are found as 
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for Rp(r) functions, and 
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(39) 

for i?o( r ) functions. Using the form of the radial wave 
functions given by Eqs. (EDI),® [Eqs. (011), (03] in 
Eqs. (|3"6 ) .([3"7 | [Eqs. ([3"5 ]) .(|3"g ]) ] we get the equation de- 
termining the energy levels of states iPfm bPool ( see ^P" 
pendix) . Solving this equation numerically we have found 
energy level positions of holes of the heavy and light hole 
types. 

The dependence of positions of the hole energy levels 
on the nanocrystal radius, their types and degeneracies 
are presented in Fig. [3] for several lowest hole levels. One 
can see that the hole level structure is more dense in com- 
parison with the electron level structure. This can lead 
to important differences in behavior of "hot" electrons 
and holes22. 



IV. COULOMB SHIFT 

The Coulomb interaction leads to a decrease of the 
recombination energy of an electron-hole pair—. This 
interaction should be considered taking into account the 
" image charge" effects appearing because of the dielectric 
constant difference at the quantum dot boundary. On the 
other side the dielectric constant mismatch leads to inter- 
action of a charged particle with its own image. Result- 
ing polarization self-energy correction increases single- 
particle energies of electrons and holes (counted upwards 
from the bottom of the bulk conduction band and down- 
wards from the top of the valence band, respectively). 
If the dielectric constant changes discontinuously at the 
quantum dot boundary the polarization self-interaction 
diverges there. This is not an immediate problem for 
models of quantum dots assuming infinite energy barri- 
ers for electrons and holes at the quantum dot boundary 
because particle wave functions vanish there and the self- 
energy correction induced by the self-interaction remains 
finit e) 24 ' 25 ' 26 . However, as a consequence of finite energy 
barriers particle wave functions are finite at the quantum 
dot boundary (see Fig. ^ and self-energy corrections be- 
come infinite. 

In order to remove these unphysical divergences one 
has to take into account that the position-dependent di- 
electric constant should change smoothly between the 
values corresponding to the quantum dot core and the 
surrounding material on the scale of the interatomic dis- 
tance. This can be done in a simple and intuitive way 
by "regularizing" directly the self-interactio n 27 ' 28 . A 
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FIG. 4: (Color online) Dependence of the ground-state 
electron-hole recombination energy as a function of diameter 
of nanocrystal (solid line). Dashed line shows the same energy 
without taking into account the exciton shift. For compari- 
son, the experimental data obtained from photoluminescence 
spectr a 32 ' 33 ' 34 ' 35 measured for Si nanocrystals inside Si02 are 
presented. 



more physical and controllable, however more demanding 
way is to solve the problem with the smooth position- 
dependent dielectric constant. Calculations of Ref. 
show that the total correction to the electron-hole pair 
energy introduced by the Coulomb and "image charge" 
interactions depends strongly on the width of the tran- 
sition region between two values of the dielectric con- 
stant. For reasonable values of the transition region 
width around the interatomic distance the corrections 
introduced by the electron-hole Coulomb interaction and 
the polarization self-energy corrections cancel each other 
to high extent. Therefore, the overall Coulomb correc- 
tion to the energy of the electron-hole pair is small and 
electron-hole recombination energy is pretty accurately 
given by the sum of the band gap with the electron 
and hole single-particle quantization energies. This con- 
clusion is also supported by GW calculations for small 
hydrogen-terminated silicon nano crystals 3 ^. 

An estimate (rather an upper limit estimate^) of the 
overall Coulomb shift Vc can be deduced using single- 
particle wave functions corresponding to infinitely high 
energy barriers. We have calculated Vc under these as- 
sumptions treating the Coulomb correction as a pertur- 
bation to the single-particle Hamiltonia n 26 ' 31 . The re- 
sult for the ground-state electron and hole can be given 
then in a simple form: Vc — — 1.54e 2 /(Ksii?nc), where 
Ksi is the dielectric constant of Si and it was taken into 
account that the surrounding material has an approx- 
imately three times smaller dielectric constant than Si 
(one should notice that Ref. [26| gives a different numeri- 
cal constant computed incorrectly by us). The excitonic 
shift calculated in the same way for higher states of elec- 
trons and hole is on the same order and smaller. 

In Fig. [4] the dependence of the exciton energy as a 
function of diameter of the nanocrystal is presented for 
the ground state. This dependence is also shown when 
the excitonic energy shift is accounted for. The theory is 
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FIG. 5: (Color online) Probabilities P r , a r of radiative transi- 
tions between the ground electron and hole states: assisted by 
emission of a TO-phonon (dashed line) , an LO-phonon (dash- 
dot line) and their sum (thick solid line) as well as probability 
of direct (zero-phonon) transition (dot line), as functions of 
nanocrystal diameter. Experimental points 3 - 5 - are shown as 
well. 



compared with the experimental data obtained from the 
photoluminescence spectr a 32 ' 33 ' 34 ' 35 . 



V. RADIATIVE RECOMBINATION 

We have produced calculation of the probabilities P r 
of radiative recombination assisted by emission of an op- 
tical transverse phonon (with energy 57.5 meV) as well 
as longitudinal one (55.3 meV). These channels of radia- 
tive transitions dominate in bulk Si. The details of the 
calculation will be published in a separate paper. The re- 
sults of calculations for the ground exciton state [P r ,gr) 
are presented in Fig. [5j The probabilities of radiative 
transitions involving excited states P r have similar de- 
pendences on nanocrystal size being of the same order of 
magnitude (e.g. for transition from the second electron 
state to the first hole state P r ks 0.8P r , gr )- In Fig. [5] 
the result of calculations of direct (zero-phonon) radia- 
tive transition for the ground exciton state is presented 
as well. Such a transition becomes possible for confined 
carriers but one can see that the oscillator strength is no- 
ticeably less for the dots with diameter larger than 2 nm. 
This is a well-known experimental fact—. One should 
notice that the shown probability of the direct radiative 
transition has been calculated as an average value over 
the nanocrystal size distribution in order to achieve the 
acceptable convergence of the numerical integration and 
to avoid strong oscillations of the resul t 37 ' 38 . One can see 
from Fig. [5] that our results reproduce the experimental 
data on radiative lifetimes in Si/Si02 nanocrystals^ 5 - very 
well. As far as raditive transitions most probably take 
place together with a phonon emission the exciton band 
gap derived from the photoluminescence spectra should 
generally be lower than the calculated one by the amount 



7 



of the phonon energy (cf. Fig. [4]) . 

VI. CONCLUSION 

We have calculated the wave functions and the energy 
levels of confined carriers in Si quantum dots inside a 
SiC>2 matrix as functions of the dot diameter. It has 
been shown that for small quantum dots [d < 2.5 nm) 
the energy spacing between neighboring electron and hole 
levels is of the order of hundreds of meV, and for elec- 
trons they are larger than for holes. Such energy spac- 
ings are also larger than the energy level splittings due 
to different mechanisms 9 , which are not accounted for 
in this paper. Therefore, the single-phonon relaxation 
between the levels becomes impossible and the time of 
relaxation of "hot" carriers to the ground state should 
increase. The calculated recombination energies of an 
electron-hole pair in the ground state and the probabili- 
ties of radiative interband transitions between the ground 
electron and hole levels are in good agreement with ex- 



perimental data. Comparison of the electron and hole 
levels calculated by our method and by state-of-the-art 
computational methods^ could lead in the future to a 
more rigorous formulation of the boundary conditions in 
the framework of the multiband effective mass approxi- 
mation. This would allow to produce rather quantitative 
calculations of processes involving electrons and holes in 
Si nanocrystals with a reasonable computational effort. 
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APPENDIX: EQUATIONS FOR QUANTIZATION 
ENERGIES OF HOLES 

Equations determining the energy quantization levels 
for holes of mixed, heavy and light holes are 
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I 



(A3) 



respectively. Here [i = B/A, v = A a /A, (3 is defined determined by Eq. (fig)) and Eq. (f23|) . respectively, 
by Eq. (|T7j) , A and k are functions of the hole energy E 
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